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Bound, strange, neutral superheavy nuclei, stable against strong decay, may exist. A model 
effective field theory calculation of the surface energy and density of such systems is carried out 
assuming vector meson couplings to conserved currents and scalar couplings fit to data where it 
exists. The non- linear relativistic mean field equations are solved assuming local density baryon 
sources. The approach is calibrated through a successful calculation of the known nuclear surface 
tension. 
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I. INTRODUCTION 



An inspection of the semiempirical mass formula 
(SEMF) reveals that the single largest limiting factor in 
the creation of very large nuclei is the Coulomb repulsion. 
One way to overcome this barrier is to include hyperons 
in nuclei 0, |, ^. Consider the hyperons A°, E(='=''^\ and 
5;(-.o)_ ^]^g lightest hyperon, the A, has a negative bind- 
ing energy in nuclear matter Q and decays weakly into 
nonstrange matter. The E's appear to have a repulsive 
nuclear potential ^, 0|. Next in mass are the S's; ex- 
perimental evidence suggests that the binding energy of 
a single S in nuclear matter is negative l^], ^, |l^ . In ad- 
dition, the reaction 2A — > -I- S becomes energetically 
favorable for some critical number of A's in the nuclear 
medium As a result, we expect that for large systems 
the addition of A's and S's is desirable, but the inclusion 
of S's would have little or no positive effect. Therefore we 
consider matter composed solely of p, n. A, 'EP , and S^. 
The inclusion of the S~'s offsets the Coulomb repulsion 
of the protons; this potentially allows for the creation of 
arbitrarily large nuclei by diminishing the importance of 
the Coulomb term in the SEMF. To minimize the effect 
of the Coulomb term we investigate this class of nuclei 
such that Q = 0. Due to the fact that these nuclei are 
stable against strong decay, they consequently decay on 
weak interaction timescales enhancing the potential for 
their detection. The purpose of this paper is to model 
the surface structure, and acquire the surface energy from 
the calculated SEMF, for this class of nuclei. 

To accomplish this, we must solve the nuclear many- 
body problem. Using quantum hadrodynamics (QHD), 
we construct an effective lagrangian density invariant un- 
der SU (2) L (E) SU{2) R symmetry using hadrons as de- 
grees of freedom Isodoublet baryon fields, t/i, 
are coupled to the meson fields 0, V, p, and tt; the me- 
son fields correspond to an isoscalar, Lorentz scalar and 
vector, and an isovector, vector and pseudoscalar, fields 
respectively In the relativistic mean field limit, the 
meson fields become classical fields and the sources be- 
come expectation values. The pions have no mean fields 
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in a spherically symmetric system. Then the lagrangian 
density is converted into a hamiltonian density; since it 
is explicitly time independent the hamiltonian density 
is equivalent to the energy density. The baryon mass 
is determined by solving the scalar field equation self- 
consistently at each point |]l2[ The coupling con- 

stants are fit to reproduce experimental values of vari- 
ous ordinary nuclei; specifically the parameter sets NLC 
and Ql, which include nonlinear scalar couplings, are 
used [|ll|, . This theory, when solved in the relativistic 
Hartree approximation, has had great success in predict- 
ing the bulk properties and shell structure in ordinary 
nuclei 

To investigate the structure of the many-particle super- 
heavy nuclei of interest here, we directly use the approach 
of density functional theory (DFT). DFT tells us that the 
exact ground-state density is acquired by minimizing the 
energy functional jlj, The effective lagrangian dis- 
cussed above provides a lowest order density functional. 
The scalar and vector fields here play the role of relativis- 
tic Kohn-Sham potentials In order to calibrate 
our approach, we calculate the ground-state densities of 
ordinary finite nuclei with N = Z. To model finite nuclei 
we must add spherically symmetric spatial variations of 
the meson fields to the lagrangian density. The source 
terms are evaluated using a local density approximation; 
at every point within the nucleus the baryons are as- 
sumed to be a local Fermi gas with states filled up to fci?. 
We acquire the scalar mean field equation by minimizing 
the energy functional with respect to the scalar field; a 
similar approach yields the vector mean field equation 

The nonlinear scalar field equation is solved as a fi- 
nite difference equation utilizing a shooting method. The 
boundary conditions are determined by assuming the 
baryon density vanishes at the surface and solving the 
linear scalar field equation outside |Q. However, these 
boundary conditions are exact only in the linear case; a 
correction term must be added to compensate for the ef- 
fects of nonlinear terms in the scalar field equation. The 
total energy is minimized with respect to the local Fermi 
wave number, while keeping fixed the baryon number, B. 
The constraint of fixed B is incorporated with a Lagrange 
multiplier, which is the chemical potential. The result- 
ing constraint equation states that the chemical poten- 
tial must be constant throughout the nucleus This 
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approach is more sophisticated than a simple Thomas- 
Fermi method because we self-consistently solve for the 
source terms at each point. With the calculated binding 
energy and baryon number for finite nuclei, we fit the first 
two terms in the SEMF, the bulk and surface energies, 
for nuclear matter. These are in good agreement with 
the known experimental values ||T2|, thereby validating 
our approach. 

Now we add in hyperons; however, existing experimen- 
tal data requires that some assumptions be made. First, 
we couple universally to the conserved baryon and isovec- 
tor currents. Second, a different scalar coupling is used 
for each baryon. The scalar coupling for the A's is fit such 
that the binding energy of a single A in nuclear matter 
is -28 MeV [||. However the binding energy of a single 
S is relatively uncertain, values appearing in the litera- 
ture range from -40 to -14 MeV. Recent experiments with 
light nuclei suggest that the value lies on the less bound 
side of this range 0; however, it may be more deeply 
bound for heavy nuclei ||l8|] . As a result, a number of 
values for the S scalar coupling are investigated. Third, 
we continue to utilize the parameter sets for ordinary nu- 
clear matter, NLC and Ql, to generate the nucleon and 
non-linear scalar couplings. 

The addition of new baryons in the theory only re- 
quires the inclusion of new source terms in the effective 
lagrangian density. We investigate a specific sector of the 
theory by imposing the restrictions Q = and |S|/B = 1, 
where S is the total strangeness. We also assume an av- 
erage cascade mass. Note that the minimum binding en- 
ergy always occurs such that there are equal numbers of 
n and p (and consequently equal numbers of "EP and S^); 
therefore the symmetry term in the SEMF is rendered ir- 
relevant. Since there is only one chemical potential, the 
reactions 2A ^ N + 'E. and n + 'Ef ^ p + S~ are in equi- 
librium. Again, using DFT to model finite nuclei, we 
now investigate the role of the first two terms in the cal- 
culated SEMF. Also, by determining the baryon density, 
we acquire the structure of the surface. 

At T = and normal nuclear densities, the mass dif- 
ference between strange and nonstrange quarks is less 
then the Fermi energy of massless nonstrange quarks. 
This opened the possibility that strange quark matter 
composed of u, d, and s quarks might be stable against 
strong decay and perhaps even absolutely stable ^ . 
These systems are characterized by small charge fraction 
Q/B ~ and large strangeness fraction |S|/B ^ 1. The 
plausibility of bound strange matter has also been ex- 
plored in the hadronic sector. Theoretical investigations 
of multiple A hypernuclei indicate that they are bound 
and stable against strong decay. These studies produced 
systems with binding energies as low as -9 MeV corre- 
sponding to |S|/B ~ 0.2 H, H, H, ||. 

Gal ct. al. suggested that discussions of matter com- 
posed of n, p, and A must also include S° and S~ due 
to the fact that the reaction 2A ^ -f S is energet- 
ically favorable for some critical number of A's in the 
nuclear medium They investigated possible configu- 



rations of bulk matter in the relativistic mean field ap- 
proach, suggesting binding energies per baryon as low 
as -25 MeV with a large strangeness fraction; finite nu- 
clear calculations were also preformed ||, g6[ |2^. In 
addition, Gal et. al. fit to a generalized SEMF using a 
Fermi gas model [^8[ |9|. Extrapolating from the ordi- 
nary SEMF they estimate the bulk and symmetry terms, 
while leaving the Coulomb term unchanged. The surface 
energy is scaled as inversely proportional to the aver- 
age baryon mass, yielding a value of 15 MeV. Stoks and 
Lee challenged these findings using a many-body theory 
with baryon-baryon potential models. These potentials 
were developed using an SU (3) extension of the Nijmegen 
soft-core potentials |3^, ^ In contrast, the lat- 

ter found that this type of matter is only slightly bound, 
E/B ~ —3 MeV or less A quark-meson coupling 

model produced a minimum binding energy of -24.4 MeV 
with |S|/B ^ 1.38 H. The effect of adding hyperons 
has also been explored in application to neutron stars 



In this paper we first construct the effective lagrangian 
of the theory following the methodology of Furnstahl, 
Serot, and Tang jll, 13, 40j. For ordinary nuclear matter, 
we model both infinite matter and finite nuclei of various 
radii using QHD. Initially, it was assumed that, because 
of the relatively large value of vector meson mass, the 
derivatives of the vector field were negligible how- 
ever, iteration on the vector field produces a significant 
effect on the calculation. Also, a correction to the bound- 
ary condition is required to account for the nonlinearities 
in the scalar field equation. From this calculation, we 
acquire the baryon density and scalar field for ordinary 
finite nuclei; this provides a picture of the size and shape 
of the surface. Then the surface energy is extracted by 
fitting to the SEMF. The approach is successfully cal- 
ibrated by comparing this result to the experimentally 
obtained values. Next this system is extended to cascade- 
nucleon matter, composed of p, n, "EP, and subject 
to the constraints Q = and |S|/B = 1. However, due 
to the uncertainty in the binding of a single cascade in 
nuclear matter, we examine a range of values for the 2 
scalar coupling. Again, by calculating finite systems, the 
bulk and surface terms are extracted from the SEMF. 
A's are then added to the cascade-nucleon matter, with 
the A scalar coupling fit to experiment. The bulk quanti- 
ties change little, however, from the cascade-nucleon case. 
An investigation of possible hyperon-hyperon interaction 
is conducted by coupling a $ meson to the conserved 
strangeness current; we allow the $ coupling to increase 
until the many-body system is no longer bound in order 
to find the maximum allowable value of this coupling. 



The organization of this paper is as follows. In section 
2, the effective lagrangian is constructed and the sys- 
tems of equations are derived. In section 3, we present 
the methodologies used to solve these systems described 
in section 2. In section 4, the results for nuclear and 
hyperon- nucleon matter are discussed. 
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II. THEORY 

Following the QHD approach of Furnstahl et. al. 
0, |l^, we construct an effective lagrangian density 
using hadronic degrees of freedom that remains invariant 
under SU{2)l ^ SU{2)r symmetry. We will use this 
lagrangian density to model both infinite and finite sys- 
tems; to start with only systems of nucleons with N = Z , 
which we refer to as nucleon matter, are considered. The 
full lagrangian density (with ?i — c = 1) is given by 



1 

2 ^ M 



/a0 



3! 



4! 
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(1) 



Here ■0 



is the isodoublet baryon field, 7^ 



{id(i,(i) are the gamma matrices, and ijl,v — 1,...4. 
The conventions of are used. The Lorentz scalar 
meson field, (j), is coupled to the scalar density ipil) and 
VJi, the Lorentz vector meson field, is coupled to the con- 
served baryon current irp^^ip; their respective coupling 
constants are gs and g^- The masses of the nucleon, 
scalar, and vector fields are denoted by M, rus, and ruy 
respectively, k and A are constants that determine the 
strength of the nonlinear scalar couplings. The field ten- 
sor is defined as 



p — 



(2) 



Notice the p-meson terms in Eq. have been sup- 
pressed. The source term contributed by the p-meson 
depends 0x1 N — which vanishes for the systems under 
consideration. We now employ relativistic mean field the- 
ory (RMFT). In this case the meson fields are replaced 
by their classical, or mean fields 



Mr) 



(3) 
(4) 



which are time- independent. For the purposes of this 
paper, we restrict ourselves to spherical symmetry. The 
source terms are now replaced by their expectation val- 
ues. In DFT, the classical fields play the role of Kohn- 
Sham potentials |l^, |l^ . Incorporating both Eq. (||) and 
RMFT, our lagrangian density becomes 
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(5) 



where the effective mass is defined as 
M* EE M - gs(l)o 
Now the hamiltonian density is given by 

dtp 



(6) 



Hix^) 



n 



at 



dx 



(7) 



where pB = ip'^fp is the baryon density and the canonical 
momentum density is 



n 



dC 



didip/dt) 



(8) 



If one assumes that a nucleus is a local Fermi gas filled 
up to some kp(r) at every point, the source terms take 
the form 



PB{r) = (ipy) 



t,/,\ - _L 



(27r)3 



(9) 



(10) 

where p = — iV and 7 is a degeneracy factor. Since the 
hamiltonian density is explicitly time- independent, it is 
equivalent to the energy density (with 7 = 4 for nucleon 
matter) 



£{r) 



K A 

+gvPBVQ + + —(j)'^ 



+ 



(27r)3 



(fc2 + M*2)l/2rf3^ 



The total energy and baryon number are 

E = / £{r)(fr; ( pB{r)(fr 



(11) 



(12) 



The energy density above provides a lowest order den- 
sity functional. DFT tells us that minimizing the exact 
energy functional yields the exact ground-state density. 

To conduct infinite nucleon matter calculations, we ne- 
glect spatial variations in the meson fields; the resulting 
energy density is 

^ = \ml(t)l + + ^(j)^ + gvPsVo - ^mlV^ 



(27r)3 



{k" + M*^Y/^d'k 



(13) 
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By minimizing the energy functional with respect to the 
scalar field, the scalar mean field equation is determined. 
The vector mean field equation is similarly derived as an 
extremum of the energy functional. These equations are 



2m? 



A 



9v 



9s 



tPb 



where the scalar density is given by 



Ps = 



(27r)3 



M* 



(fc2 + Af*2)l/2 



(14) 
(15) 

(16) 



The solution to these equations is discussed in section 3. 

We now turn our attention to finite nucleon systems. 
We retain the spherically symmetric spatial variations 
in the meson fields and therefore require the full energy 
density in Eq. (|lTl). The meson field equations, acquired 
in the same manner as above, are 



(V2 - ml)Vo 



A 
6' 



-9sPs 



-9vPb 



(17) 
(18) 



Note that for spherically symmetric systems, the lapla- 
cian becomes 



2 92 2d 

V = 7T H 

Using a Green's function, the solution to Eq. ( p^ ) is 

-tUtj \x — y\ 



(19) 



9vVo 



f / d^VPBiyf-r^ 
At: J \x- y\ 

9v 



xrriv 



ydypsiy) sinh(m„x<)e" 



(20) 



(21) 



where the angular dependence has now been integrated 
out. Since the contribution of the laplacian in Eq. ( p^ is 
small compared with that of the vector meson mass, to 
a first approximation it can be neg lected @ g, 

9vVo = (22) 

However, the omitted term produces a small, but im- 
portant, contribution to the vector field; this can have a 
significant effect on the total energy and baryon number. 
It is then convenient to express the vector field as 



gvVo = ^PB + SWo 



(23) 



with SWq = gySVo- Substituting this in Eq. ( pi] ) and 
rearranging, one obtains an explicit expression for SWq 
in terms of the baryon density 



2 f -m^\x-y\ 

SWoir) = / d'ypBiy) 

in J \x- y\ 



^PB{r) (24) 



Minimization of the total energy with respect to the 
local Fermi wave number now yields the ground state of 
the system. A Lagrange multiplier is used to incorporate 
the constraint of fixed B such that 



(25) 



Since the variations of the energy density with respect to 
both the scalar and vector field vanish, they can be held 
constant in the variation of kp , the result of which is the 
constraint equation 



p^^9vVo{r) + [ep{r) + M*\r)Y/^ 



(26) 



where the Lagrange multiplier, /z, is the chemical poten- 
tial and is constant throughout the nucleus. 

On the surface r = rg, the baryon density vanishes. 
The constraint equation at the surface then yields the 
first boundary condition 



(27) 



where Eq. ( |23D has now been employed. 

To determine the second boundary condition, consider 
the solution to the linear homogeneous scalar field equa- 
tion 



-ij--ra)m.s 



(28) 



The constant gs<j)c = M — p + SWo{ro) is determined 
at the surface. Differentiating Eq. ( psf ) with respect to r 
and then evaluating at the surface, we acquire the second 
boundary condition 



dM*{r) 
dr 



= [M + <5VKo(ro)-Ai] 



1 4- msTQ 



(1 



(29) 

The solution to the scalar field equation in Eq. (2^) 
no longer holds when the nonlinear terms arc included; 
therefore, a small correction e has been included to com- 
pensate. In the nonlinear case, the scalar field equation 
is integrated outward from tq , and e in Eq. ( |29| ) is varied 
until the solution vanishes for large r. 

The calculated binding energy and baryon number of 
a series of nuclei with different radii can be fit with a 
SEMF of the form 



E 

— = fli 
B 



a2 



(30) 



where only the bulk and surface terms have been re- 
tained. The bulk constant, oi, is determined by the 
binding energy of infinite nuclear matter. Then, after 
calculating a number of finite nuclei, the surface energy, 
02, can be obtained by plotting the calculated energies 
per baryon against the calculated values of i?^^/^. The 
results of numerical methods for these finite systems are 
discussed in later sections. 
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We now extend our theory to consider systems of nucle- 
ons and hyperons. First we investigate matter composed 
of n, p, S", and S^, subject to the conditions 



g-0 

\S\/B = 1 



(31) 
(32) 



where Q and S are the total charge and strangeness re- 
spectively. These systems shall be subsequently referred 
to as cascade-nucleon {'B.N) matter. Equation ( |3l] ) re- 
stricts the system to equal numbers of p and S^; sim- 
ilarly, Eq. (32) forces the numbers of n and S~ to be 
equal. Therefore the system is now characterized by two 
Fermi wave numbers, kpp and kpn- For simplicity we 
employ an average cascade mass Ms = (M^o -I- M^-)/2. 
Since the energy density is now symmetric under the in- 
terchange of kpp and kpn, the minimum binding energy 
always occurs such that kpp — kpn. As a result, we 
can further restrict this system to a single Fermi wave 
number, kp. It is a consequence of these arguments that 
equilibrium is imposed upon the reaction 



(33) 



and the system is described by only one chemical poten- 
tial. Again we mention that the p-mesons do not con- 
tribute here for similar reasons to the nucleon case. 

Next, we must make some assumptions about the cas- 
cade couplings. Since the baryon current is conserved, 
the vector coupling is taken to be universal, = g^s- 
However, an independent scalar coupling for the cas- 
cades, gsE, is assumed. 

Consider the case of infinite BN matter. The addition 
of hyperons to the theory requires only the addition of 
new source terms. A source term of the form 



S£ 



7 



{2n) 



{k^ + AG^)'^^d^k 



(34) 



is added to the energy density in Eq. (13) where 

Mi = Me - gsE^o (35) 
In addition, a new term is included in the baryon density 



1 



(2^)^ 7o 



(Fk 



(36) 



Here ^ — A for (S°,S^) with spin up and down. Except 
for the additional source terms, the meson field equations 
remain unchanged. The new term added to the source in 
the scalar field Eq. ( [T^ is 



5ps 



As 
(2^ 



d^k 



Mi 



(P + M|2)l/2 



(37) 



where 



Equation (p6)is incorporated into the source in the vec- 
tor field Eq. (15|). The solution to these equations is again 
discussed in section 3. 

We now examine the case of finite SiV matter. The 
source terms in Eqs. (|^ and ( ^ ) are incorporated into 
t he energy density in Eq. ([III). Next the terms in Eqs. 

are added to meson field equations in Eqs. 
18D respectively. Then a new constraint equa- 



(37) and (3' 
and (T 

tion is produced in the same manner as before 

= 9vVo + \[{kl + A/*2)i/2 + (4 + M|2)i/2] (39) 

Similarly, the boundary conditions are now 
1 



M*(ro) = 



1 



-[2^ - 2(5M/o(ro) - Me + sM] (40) 



dM*{r) 



dr 



-[25Wo{ro) - 2^1 + Me 



+M]['-±^']{1 



(41) 



Consider again the SEMF in Eg (pO[) . The conditions 
imposed on EN matter in Eqs. ([si]) and ( ^2|) now jus- 
tify the elimination of the Coulomb and symmetry terms. 
Then oi is taken to be the binding energy of infinite EN 
matter; next, proceeding as before, the surface energy, 
02, can be extracted. 

Finally, we investigate a class of matter in which A's 
are added to the EN matter described above. These sys- 
tems are referred to as lambda- cascade-nucleon [KEN) 
matter. The previous restrictions do not relate the num- 
ber of A's to the number of TV's and S's; therefore a 
second Fermi wave number, fc^^A, is needed. Again the 
vector coupling is taken to be universal and an indepen- 
dent scalar coupling, gsA, is employed. Now equilibrium 
is imposed on the reactions 



n + S° ^ A + A 



A + A 



(42) 
(43) 



as well as on Eq. (p3|). As before, the system is then 
characterized by a single chemical potential. The source 
terms required for the inclusion of A's in both infinite 
and finite EN matter are 



6E = 



(27r)3 



(P + Miy/^d^k (44) 



(27r)3 



Sps = ,^ / d^k 



d^k 



Ml 



(27r)3 7o (fc2 + Mf)i/2 



(45) 
(46) 



where 



s = gsE/gs 



(38) 



MX = Ma - 9sA<l>o 



(47) 
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Number of Iterations 



FIG. 1: Convergence of the baryon number and binding en- 
ergy per baryon (boxes and circles respectively) after 9 iter- 
ations for an ordinary nucleus of ro = 15/ms, N = Z, and 
using the L2 parameter set. 







2 


ms 


m„ 




A 


QsA/gs 


L2 


109.63 


190.43 


520 


783 








.886601 


NLC 


95.11 


148.93 


500.8 


783 


5000 


-200 


.881898 


Ql 


103.67 


164.70 


504.57 


782 


4577.6 


-197.70 


.884029 



TABLE I: Parameter sets taken from |l^, nis, m^, k, and 
A are in MeV. Qsa/Qs is fit to reproduce the binding energy 
of a single A in nuclear matter 



9ss/gs 


S£/5pB - Ms 


1 


-68.9 


.95 


-51.6 


.9 


-34.3 



TABLE II: Values of the binding energy (in MeV) of a single 
H in nuclear matter for various H coupling ratios and the 
parameter set NLC. 



and 



t = 9sa/9s 



(48) 



In the case of finite A.^A^ matter, there are now two con- 
straint equations 



(50) 



The A density begins interior to the surface ro; this allows 
the 5iV boundary conditions to be used in the AS7V case. 



III. METHODOLOGY 

In this section we develop a methodology for solving 
the systems of equations discussed in section 2. In the 



case of nucleon matter the parameters gs, gv, i^is-, itt-v, M, 
K, and A must first be specified. The vector meson mass 
is defined to be the mass of the w-meson and M = 939 
MeV. The remaining constants are given by the three 
coupling sets in Table || to determine these sets the the- 
ory was fit to reproduce various properties of ordinary 
nuclear matter [|l^ The simplest set, L2, includes 
only linear terms in the scalar field. The sets NLC and 
Ql both expand the theory to include nonlinear terms. 
As a result, these sets must be fit to more properties of 
nuclear matter than L2. 

To extend the theory to systems of nucleons and hy- 
perons, specification of the constants gsA, gss, <?i,a, and 
gyS is also required. Since the vector meson is coupled 
to the conserved baryon current, we assume a universal 
vector coupling, = gy\ = g^-^. The scalar couplings, 
on the other hand, are adjusted to reproduce the bind- 
ing energies of single hyperons in nuclear matter. For 
instance, the A scalar coupling is designed to reproduce 
the binding energy of a single A in nuclear matter, ex- 
perimentally determined to be —28 MeV The values 
of gshl gs are also shown in Table ^ Unfortunately data 
on the binding energy of a single S in nuclear matter is 
uncertain. Therefore, a range of S scalar couplings is in- 
vestigated; the values used are gss/gs ~ 1, -95, and .9. 
These values correspond to the binding energies listed in 
Table 

Consider the case of infinite nucleon matter. To ob- 
tain the solution to this system, first one must specify 
kp- Both Eqs. and (|l|) are now solved for their re- 
spective meson fields. Then, using the meson fields and 
kp, one calculates the energy density in Eq. (|l^). This is 
in turn used to evaluate the binding energy per baryon, 
BE(A;f) = £/pB - M. In RMFT the medium saturates 
and BE(fci?) has a minimum; this equilibrium value, BEq, 
serves as the bulk term in the SEME, ai. This procedure 
is also applicable to infinite HA'' and AEN systems pro- 
vided the appropriate source terms are included. 

Now we discuss the methodology used for all finite sys- 
tems. First the scalar field Eq. (|l^) is converted into a 
pair of coupled first-order finite difference equations for 
[(j)o{r), (j)Q{r)]; these equations are solved using a shoot- 
ing method. To accomplish this, one fixes /i and ro; now 
the boundary conditions are uniquely determined. Since 
kpifo) — 0, we can solve the finite difference equations 
for 00 (r) and 0o(r) one step in. These solutions are sub- 
stituted back into the constraint equation from which 
kp{r), one step in, is determined. In this manner we it- 
erate in from the surface, evaluating /Os(r) and £{r) at 
every point. This process is repeated until ps and 0o be- 
come constant as r approaches the origin (or equivalently 
p'b = 't'o = ^ — 0)j achieve this by adjusting the 
chemical potential. We mention that better convergence 
is obtained by decreasing the step size. Note that ini- 
tially, e and SWq are ignored. 

Now we incorporate the two correction terms. In order 
to discuss the role of e, let us examine the second bound- 
ary condition. When nonlinear terms are introduced, the 
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FIG. 2: Binding energy per baryon for infinite nuclear matter 
with N = Z as a function of Fermi wave number. These results 
are for the coupling sets NLC (solid line) and Ql (dashed 
line). 



entire process is repeated again. After successive itera- 
tions on the vector field, B and BE = E/B — M both 
converge to their fuh solutions. This convergence is il- 
lustrated in Fig. for ordinary nuclcon matter; here, the 
parameter set L2 was used, a radius of tq — 15toq ^ was 
assumed, and 9 iterations on the vector field have been 
carried out. Similar convergence was found in all cases 
studied here. 

Finally, we consider the SEMF in Eq. (^0|) where a\ 
is defined as the BEq of infinite matter. The calculated 
values of BE and B"-'^/^ of the finite systems are plot- 
ted against each other for nuclei of various radii. Then 
using this SEMF as a linear fit, the surface energy, 02, 
is determined. The above approach to finite systems is 
first calibrated by the nucleon matter case; then it is ex- 
tended to investigate the EN and AS7V systems detailed 
in section 2. 



IV. RESULTS AND DISCUSSION 





kp 


M7M 


£/pB ' M 


L2 


1.301 


.5409 


-15.758 


NLC 


1.301 


.6313 


-15.768 


Ql 


1.299 


.5975 


-16.099 



TABLE III: Calculated equilibrium values of the Fermi mo- 
mentum (in fm~^), effective mass, and the BEq (in MeV) for 
infinite nucleon matter are shown using the couplin g se ts in 
Table ^ These numbers reproduce the results in p^]4C|] . 







kp 


M7M 


£/pB - Ma 


NLC 


1 


1.363 


.1391 


-41.343 


.95 


1.351 


.1686 


-22.747 


.9 


1.326 


.2232 


-4.741 


Ql 


1 


1.320 


.1372 


-42.728 


.95 


1.310 


.1643 


-24.013 


.9 


1.290 


.2104 


-5.812 



TABLE IV: Calculated equilibrium values of the Fermi mo- 
mentum (in fm~^), effective mass, and the BEo (in MeV) for 
infinite HV matter are shown using the coupling sets NLC 
and Ql in Table | and a range of values for gss/gs- 



small parameter, e, is included to compensate. Iterating 
the finite difference equations out from the surface, e in 
Eq. ( p9| ) is adjusted such that the scalar field vanishes 
for large r. The newly corrected boundary condition is 
then used to resolve the finite system by integrating in 
as described above. 

Next, the correction term 6W0 is added to the vector 
field. Initially Eq. (^2|) was employed; however, this is 
accurate only in the limit of a large m.^,. A small, but 
important, contribution to the vector field was omitted; 
therefore, the term SWq, defined by Eq. (24) and cal- 
culated from the previous pB(r), is included. Then the 



In this section, we discuss the application of the above 
methodologies. First we consider the results of our infi- 
nite matter calculations, starting with nucleon systems. 
Table III shows the equilibrium values kp, M*/M, and 
BEq obtained for infinite nucleon matter using the L2, 
NLC, and Ql parameter sets. The values calculated here 
reproduce those in [^3[ Then BE(fci?) is plotted for 
both the NLC and Ql sets in Fig. ||. The minimum, BEq, 
in Fig. 1^ is taken to be ai in the SEMF for each coupling 
set. This is in good agreement with the empirical bulk 
term @. 

We now turn our attention to infinite systems of nu- 
cleons and hyperons; these investigations are conducted 
using only the nonlinear cases, NLC and Ql. We begin by 
investigating infinite SV matter for the range of S scalar 
couplings mentioned above. Since the lowest mass state 
of separated baryons under the conditions Eqs. ( |3l| ) and 
(^) consists entirely of A's, BE{kp) for infinite cascade- 
nucleon matter is defined by BE(fci?) = 8/ps~ Ma- The 
calculated equilibrium values for this type of matter are 
given in Table [V. Notice that |BEo| decreases as the 5 
coupling grows weaker while the equilibrium kp remains 
fairly constant. Graphs of BFi{kp) for the NLC set and 
each S scalar coupling ratio shown in Fig. ^ illustrate this 
point. Although the equilibrium kp is roughly the same 
here as in the nucleon case, these systems contain twice 
as many baryons; as a result, the baryon density is much 
higher than in infinite nucleon matter. Also, the effective 
mass is considerably smaller, on the order of a third the 
value of the nucleon case. Again ai in the SEMF for each 
coupling ratio is taken to be the minimum (BEq) of the 
corresponding curve in Fig. |^. We also mention that the 
value of gss/gs — -886 is the lowest S coupling ratio for 
which infinite SV matter was still bound. 

Next, we consider infinite ASV systems using the same 
range of S scalar couplings quoted above. Note that as 
in the EN case, BE{kp) = £ / ps — Ma. This investiga- 




kp(fm'') 



FIG. 3: Binding energies per baryon for infinite cascade- 
nucleon matter computed relative to isolated lambdas (the 
lowest energy free baryon state for |S|/B = 1) as a function 
of the Fermi wave number using NLC. Note the left hand in- 
tercept is (Ms + Mn)/2 - Ma. The solid, long dashed, and 
short dashed lines correspond to Qss/ga = 1, 0.95, and 0.9 
respectively. 





gss/gs 


kp 


kpA 


M7M 


£/pB - Ma 


NLC 


1 


1.343 


.8665 


.1202 


-42.229 


.95 


1.319 


1.026 


.1321 


-24.704 


.9 


1.288 


1.159 


.1511 


-8.100 


Ql 


1 


1.302 


.8141 


.1213 


-43.457 


.95 


1.278 


.9831 


.1325 


-25.792 


.9 


1.247 


1.124 


.1495 


-9.095 



TABLE V: Calculated equilibrium values of the Fermi mo- 
menta (in fm~^), effective mass, and the BEo (in MeV) for 
infinite AHA'' matter are shown using the coupling sets NLC 
and Ql in Table | and a range of values for gas/gs- 
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FIG. 4: The baryon density ns(r) = pfl(r)/M^ (solid line) 
and effective mass M*(r)/M (dashed line) vs. r (in units of 
mj^) for an ordinary finite nucleus with N = Z, B = 188.869 
ro = 20/m3, and using the NLC parameter set. 





ro 




B 


E/B-M 


L2 


15 


4.7055924 


54.568 


-8.7527 




20 


4.69755554 


160.72 


-10.939 


NLC 


15 


4.6968452 


66.438 


-11.297 




20 


4.69165108 


188.87 


-12.719 




25 


4.68905238 


408.43 


-13.464 


Ql 


15 


4.6965099 


62.581 


-11.262 




20 


4.69083149 


179.76 


-12.808 




25 


4.68800689 


390.77 


-13.615 



TABLE VI: Results of finite nucleon matter for the L2, NLC, 
and Ql parameter sets and various radii. Calculations with 
L2 used 9 iterations on the vector field while 5 were used with 
NLC and Ql. The radii are in m~^, the chemical potential is 
in fm~\ and E/B - M is in MeV. 



tion produced the equilibrium values shown in Table 
Here the equiUbrium values of kp, M*/M, and BEq differ 
little from the SiV results for large S coupling; however, 
the difference becomes more pronounced as the 5 cou- 
pling decreases. In our formulation, a second Fermi wave 
number, fc^A, was included for the A's; as one might ex- 
pect, kp\ grows, and consequently the proportion of A's 
increases, as the gap between g^s and gsA narrows. The 
smallest value of the S scalar coupling for which infinite 
AS7V matter was still bound was Qss/gs = -875. 

Now we examine the results of the finite matter in- 
vestigation. To begin with, we consider the finite nu- 
cleon matter system. The calculated values of /i, B, 
BE = E/B — M as a function of rg for this type of matter 
using the L2, NLC, and Ql sets are shown in Table 0. 
As stated above, 9 iterations on the vector field were con- 
ducted on nuclei calculated using the L2 set; this demon- 
strated the convergence of the system. Subsequent finite 
nucleon matter results were obtained using 5 iterations 
which gave results for B and BE to better than 1%. The 



radii used here, tq = 15, 20, and 25 in units of mj^, 
include nuclei spanning a range of B ~ 50 — 400. As an 
example, psif) and (f)o{r) for a nucleus with tq = 20mJ"'^ 
calculated with the NLC set are displayed in Fig. ^. The 
interior of the nucleus is roughly constant in both psif) 
and (po (r) . Then the effective mass increases to near unity 
and the baryon density drops to zero at ro; this is a typ- 
ical example of the surface structure for finite nucleon 
systems. Next using the NLC set the calculated values 
of BE are plotted vs. B~^^^ in Fig. ^ for nuclei of various 
radii. Notice that the infinite matter value, BEo, has also 
been included. A SEMF of the form Eq. (|o|) is used as 
a linear fit in Fig. ||; the slope of this fit is the surface 
energy, in this case a2 = 18.0 MeV. The surfac e en ergies 
for the various coupling sets are given in Table VII along 
with the experimentally determined value ||l^ ; the values 
of 02 for both NLC and Ql show good agreement with 
experiment. The agreement between the values calculated 
with the more realistic interactions and the empirical re- 
sult for the surface energy of ordinary nucleon matter 
using this effective lagrangian and density functional ap- 
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0.2 
3-1/3 



FIG. 5: Fit to the calculated SEMF for ordinary nuclear mat- 
ter with N = Z and the NLC couplings. The surface energy 
is given by the slope of the curve, here a2 = 18.0 MeV. 





a2 


L2 


26.510 


NLC 


18.008 


Qi 


19.107 


Expt. 


17.8 



TABLE VIL Calculated values of the surface energy (in MeV) 
for nucleon matter using the parameter sets in Table |^ The 
experimental value is also included ]l^ . 



proach gives us some confidence in our exploratory study 
of the surface structure of strange superheavy nuclei. 

Now we investigate finite EN matter for the values of 
gss/ds quoted above. Since the best fit to both infinite 
and finite nucleon matter was obtained with NLC, we 
use this set exclusively in the following discussion. The 
values of /i, B, and BE = E/ B — M a obtained for finite 



SiV matter are given in Table VIII. For the same reasons 
as in the infinite SiV case, the binding energy per baryon 
is redefined as BE = E/B — Ma. Note that due to a 
slower rate of convergence, these systems were calculated 
using 9 iterations on the vector field. Also the radii, 
ro = 10 and 15 in units of mj^, were used here; this 
includes nuclei with baryon numbers ranging from B ^ 
30—200 depending on the S coupling. It is also important 
to mention that for the coupling ratio gss/gs = -9, the 
nuclei were unbound for our choice of radii. For one 
nucleus of this type, Fig. || shows the plots of both psir) 
and (/)o(^); this is for a nucleus with gss/gs = 1 and tq = 
15mJ^. Notice that the baryon densities in the interior of 
the nucleus are much larger than those in nucleon matter; 
also, the effective mass drops to less than a third of the 
nucleon matter value in the interior. The result is a much 
higher total B for a fixed rp. Another feature of note 
is the surface structure; here the width of the surface 
has decreased relative to the previous case. For each of 
the S scalar coupling ratios, the calculated values of BE 



are plotted vs. B~^^^; these plots are overlaid in Fig. 
0. The infinite matter values are also included. As in 
nucleon matter, the SEMF is used as a linear fit, one for 
each 5 coupling ratio, from which the surface energy is 
determined. The values of a2 are given in Table 



IX 



As mentioned in section 2, precise data on the bind- 
ing energy for a single S in nuclear matter is unavailable. 
The values appearing in the literature range from -14 to 
-40 MeV. This necessitated that a number of S scalar 
couplings be investigated. Now that the surface energies 
have been acquired for the S couplings, they are plotted 
vs. these coupling ratios in Fig. g. A hnear interpolation 
is used between the points; this is extended by extrap- 
olation into the region which corresponds to values of 
the binding energy of a single S appearing in the litera- 
ture. We feel confident in the neglect of the A's over this 
region because preliminary investigations of finite A2iV 
matter show that the BE and B change little from EN 
matter. In Fig. ^ the baryon density of a ASiV nucleus 
of ro ~ lOmJ^ is shown; notice that the A density begins 
interior to the surface and is comparatively much smaller. 

A preliminary calculation was also conducted with a 
$ meson coupled to the conserved strangeness current. 
This simulates a repulsion between like strange particles. 
In order to test the size of this interaction which could be 
tolerated, the $ coupling was increased until the system 
was no longer bound. The values for which this occurred 
for infinite EN matter are g<s>/gp — .68387, .50902, and 
.23525 corresponding to gss/gs — 1, -95, and .9 respec- 
tively. 

It should be mentioned that the conditions Q = and 
|S|/B = 1 were introduced to eliminate the coulomb and 
symmetry terms from the SEMF. Since both charge and 
strangeness are conserved quantities in the strong inter- 
action, these conditions are unaffected by strong (and 
electromagnetic) reactions in the system. As a result, this 
is an intrinsically interesting case and the calculation is 
simplified by the need for only a single Fermi wave num- 
ber. Of course, experimental processes could produce an 
arbitrary Q and |S|/B. Therefore it is of interest to es- 
timate how much our results might be modified as these 
conditions are relaxed. 

The SEMF has been generalized to include both nu- 
cleons and hyperons in [ p6| , |2^. The generalized 
SEMF proposed by Dover and Gal contains additional 
contributions to the bulk and symmetry energies [ p6[ . 
Their SEMF, with their parameter set 1, can be used 
to estimate how much the quantity E/B changes as one 
moves away from the conditions Q = and |S|/B — 1. 
The additional terms result in |(5E/B| < 5 MeV for the 
range 1/2 < |S|/B < 5/4 and arbitrary Q. If one makes 
the rough assumption that the calculated energy could 
change by this amount, the surface energy extracted from 
plotting E/B vs. B~^/^ could change by up to 30% (this 
is undoubtedly an overestimate) . Calculations with arbi- 
trary Q and |S|/B are more difficult. Work is in progress 
to examine some of these systems. 

It is also of some interest to consider the experimental 



manifestations of these nuclei. A number of experimental 
searches for strange matter have been conducted, exam- 
ples of which are |41, 42, ^ all have yielded neg- 
ative results. Characteristically, the systems considered 
here are stable against strong decay but unstable against 
weak decay. Therefore, their lifetimes are on the order of 
the weak interaction timescale, or 10"^*^ s. They will 
experience strangeness changing weak decays, such as the 
decay modes A -I- N ^ N -t- N and S + N ^ A + N. We 
expect heavy ion collisions or supernovae to be possible 
production sources for these systems. Of course, multi- 
strange baryon systems would have to make transitions 
to the ground state for the present calculations to be ap- 
plicable. The actual rate of production of nuclei of the 
type considered here is a question that goes well beyond 
the scope of the present paper. 

One physical consequence of our results is that the 
minimum size of these objects is predicted. This value 
is obtained by setting E/B = in the SEMF and then 
solving for B. The minimum baryon numbers derived 
from our calculations are 5.4 and 14.6 for gsE/gs = 1 
and .95 respectively. However, shell structure becomes 
more important in the region of small E/B. As a result, 
Hartree calculations are more reliable here. Research is 
in progress to more accurately estimate this number. 

In conclusion, we have used this lagrangian density and 
DFT approach to model nucleon, EN, and AEN matter. 
First, this approach was used to conduct infinite matter 
calculations. In the case of infinite nucleon matter, we 
reproduced the equilibrium values from previous works 
0. The theory was then extended to infinite :^N 
and AEN systems; the equilibrium values were also found 
for these types of matter for various E scalar couplings. 
The BEq in each case was taken to be the bulk term in 
the SEMF. Then we turned our attention to finite sys- 
tems. For the nucleon matter case, nuclei of various radii 
were calculated. Plotting BE vs. B~^/^ and using the 
SEMF as a linear fit, the surface energy was extracted. 
The values for NLC and Ql were in good agreement with 
the experimental value |l^. This gave us the confidence 
to extend the finite theory to EN systems. Calculations 
of psir) and 0o('') gave us a description of the surface 
structure of these strange superheavy nuclei. Again plot- 
ting BE vs. B-i/3 and using the SEMF as a linear fit, 
the surface energy was extracted for a given E coupling 
ratio. Then using a linear extrapolation in a plot of a2 vs. 
9sE/gs, one can ascertain the surface energy for coupHng 
ratios which correspond to the single E binding energies 
in the literature. 

I would like to thank Dr. J. D. Walecka for his sup- 
port, advice, and for independently verifying much of the 
numerical work in this paper. This work was supported 
in part by DOE grant DE-FG02-97ER41023. 
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FIG. 6: Baryon density ns(r) = pB{r)/M^ (solid line) 
and effective mass M*(r)/M (dashed line) vs. r (in units of 
mj^) for a nucleus composed of nucleons and cascades with 
ro = 15/ms, B — 164.918, and gss/gs ~ 1 subject to the con- 
straints Q = and |S|/B — 1. These results were obtained 
using the NLC parameter set. 





9ss/gs 






B 


E/B - Ma 


NLC 


1 


10 


1.1604726 


48.459 


-21.416 






15 


1.154255303 


204.64 


-28.944 




.95 


10 


1.178092 


32.866 


-5.4176 






15 


1.17206095 


164.92 


-12.241 



TABLE VIII: Results for finite HA matter for the NLC pa- 
rameter set and a number of radii. These calculations used 
9 iterations on the vector field. The radii are in m~^, the 
chemical potential is in units of M, and E/B — Ma is in MeV. 




FIG. 7: Binding energy vs. B"^'''^ for matter composed of 
equal numbers of cascades and nucleons for the NLC coupling 
set. The upper and lower curves correspond to gss/gs ~ .95 
and 1 respectively. The surface energy is just the slope of 
these lines. 
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FIG. 8: Linear fit to tlie surface energy, a2, vs. scalar cascade 
coupling ratio (pss/ffs) for cascade- nucleon matter assuming 
Q = 0, |S|/B — 1, and neglecting A's. 





gss/gs 


02 


NLC 


1 


72.687 




.95 


55.643 



TABLE IX: Values of the surface energy (in MeV) for EN 
matter using the NLC parameter set from Table ^ 
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FIG. 9: Baryon densities for a finite system of nucleons, cas- 
cades, and lambdas with B — 146.248 for the NLC coupling 
set and Qse/Qs = 1. The total baryon density, the total den- 
sity of cascades and nucleons, and the lambda density are 
shown by the solid, long dashed, and short dashed curves 
respectively. Notice that the lambda density is finite only 
interior to the surface. 
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